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Abstract 

We investigate the effects of delayed Raman response on pulse dynamics in massive multichannel 
optical fiber communication systems. Taking into account the stochastic nature of pulse sequences 
in different frequency channels and the Raman induced energy exchange in pulse collisions we show 
that the pulse parameters exhibit intermittent dynamic behavior, and that the pulse amplitudes 
exhibit relatively strong and long-range correlations. Moreover, we find that the Raman-induced 
cross frequency shift is the main intermittency-related mechanism leading to bit pattern deterio- 
ration and evaluate the bit-error-rate of the system. 

PACS numbers: 42.81.Dp,42.65.Dr,42.81.-i,05.40.-a 
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The dynamic evolution of coherent patterns in the presence of noise and nonhnearities 
is a rich and complex subject that is of major importance in many areas of physics. Fiber 
optics communication systems, which employ optical pulses to represent bits of information, 
serve as an excellent example for systems where noise and nonlinear effects can have an 
important role in the dynamics of coherent patterns 1]. It is known that the dynamics of 
the parameters characterizingthe pulses in optical fiber transmission systems can exhibit 
non-Gaussian statistics 

mm 

. Yet, it is commonly believed that the statistics of the 
optical pulse parameters is very different from the intermittent statistics encountered in 
strong nonlinear phenomena such as turbulence and chaotic flow. (For a review of inter- 

n 

mittency in the context of turbulent flow, see Ref. [6|). In this Letter we present results 
that contrast this common belief and show that the parameters of optical pulses can exhibit 
intermittent dynamic behavior in massive multichannel transmission systems. Furthermore, 
we demonstrate that this intermittent dynamics can have important practical consequences 
by leading to relatively large values of the bit-error-rate (BER) characterizing the system 
performance. 

We consider conventional optical solitons as an example for the pulses carrying the infor- 
mation and focus attention on the effects of delayed Raman response on the propagation. 
The main effect of delayed Raman response on single-soliton propagation is the self fre- 
quency shift, which is due to energy transfer from higher frequency components of the pulse 
to its lower frequency components [3, 1^ . The main effect of a single two-soliton collision in 
the presence of delayed Raman response is an energy exchange between the collidin g p ulses 
(Raman induced cross talk), which leads to a change in their amplitudes 
In addition, the frequencies of the two solitons also change as a result of the collision (Raman 



induced cross frequency shift) 
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The combined effect of Raman scattering and randomness of soliton sequences in multi 



channel transmission systems was considered in Refs. \l4. Il5l|. where it was found that the 
soliton amplitude has a lognormal distribution. It was also shown that the distribution of 
the cross frequency shift in a two-channel system is lognormal and that the first two normal- 
ized moments of the self frequency shift grow exponentially with increasing distance. Even 
though these studies implied intermittent dynamic behavior for the soliton amplitude, it was 
not clear whether the other soliton parameters exhibit similar dynamic behavior in a general 
multichannel system. Moreover, the effect of the amplitude intermittent behavior on the 



2 



main measure of system performance, the BER, was not addressed at all. In this Letter we 
study in detail the intermittent character of soliton dynamics. We show that the normalized 
nth moments of the self and cross frequency shifts are exponentially increasing with both 
propagation distance z and v? ^ i.e., the self and cross frequency shifts exhibit intermittent 
dynamic behavior. We find that the nth order equal- 2; amplitude correlation functions have 
similar dependence on z and v? with a typical correlation time which is much larger than 
the time slot width, thus showing that the system exhibits relatively strong and long-range 
correlations. Furthermore, we find that the cross-frequency shift is the main intermittency- 
related mechanism leading to an increase of the BER, and calculate the 2;-dependence of the 
BER for different channels. 

Propagation of pulses of light through an optical fiber in the presence of delayed Raman 
response is described by the following perturbed nonlinear Schrodinger equation [jj: 

i(9^^ + (9i2^ + 2|^|2^ = -ejj^9j|^|l (1) 

In Eq. is the envelope of the electric field, z is the position along the fiber, t is time 

in the retarded reference frame and the term —eB^dt\^\^ accounts for the effect of delayed 
Raman response . When 6/2 = 0, the single-soliton solution of Eq. in a given frequency 
channel /3 is described by ^/3(t, z) = ■?7^exp(zx^) cosh~"'^(a;^), where X/s = rjf^ {t — y/^ — 2/5z) 
and X/3 = + Pi^ — Up) + {rj^ — z, and ap, rfp and yp are the soliton phase, amplitude, 
and position, respectively. 

Consider a single collision between a soliton in the reference channel (/3 = 0) and a soliton 
in the [3 channel. We assume that |/?| ^ 1, ^ 1 and ^ 1/1/5 1, which is the typical 
situation in current multichannel transmission systems. Focusing attention on changes in the 
parameters of the reference channel soliton, for example, one finds that the most important 
effect of the collision is an 0(e/j) change in the soliton amplitude Il^: 

A?7o = 277o?7/3Sgn(^)eR. (2) 



The main effect of the collision in order en/ (3 is a frequency shift given by = 
~(8'7o^/3^r)/(3|/9|). Since e^^ <^ Vl/^l ^ "^^ neglect effects of order e|j and higher. 

We now describe propagation of a soliton in the reference channel undergoing many 
collisions with solitons from all other frequency channels in a system with 2N + 1 channels. 
We assume that the amplitudes of the latter solitons are all equal to 1. The stochastic 



character of soliton sequences in different channels is taken into account by defining discrete 
random variables Qj, which describe the occupation state of the jth time slot in the ith 
channel: Qj = 1 if the slot is occupied and otherwise. Therefore, the nth moment of Qj 
satisfies (Q) = s, where s is the average fraction of occupied time slots, assumed to be the 
same in all channels. We also assume that the occupation states of different time slots are 
uncorrelated: {Qjd'j') = ii i ^ i' and j ^ j'. We denote by A/? the frequency difference 
between neighboring channels and by T the time slot width. We assume that the change in 
?7o over the interval Azc^^ = T / (2 A/5), traveled by the reference channel soliton while passing 
two successive time slots in the nearby channels, is small. Using Eq. (jSj) and summing over 
all collisions occurring at Az^"* we arrive at 
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Az^ 



(1) 



''?ar^ E^g^(A)Ec.> (3) 



Z=Zk j_^Q 



where k — 1 and k are the indexes of the two successive time slots in the i = —1 channel, 
Zfc — -Zfc-i = Azi^\ and the outside sum is from — to A^. We decompose the disorder into 
an average part and a fluctuating part: Qj = s+Qj, where {Qj) = 0, {(ijQ'j') = s{l — s)6ii'Sjj', 
and 6ii' is the Kronecker delta function. Substituting Qj = s + Qj into Eq. © and going to 
the continuum limit we obtain 



1 drio 4:seRAp 



^sgn(A)N+2e«e^°)(^;iV), (4) 



Vo dz T 
where the continuous disorder field ^^^\z; N) is 

ki 

^) = -71) E E (5) 

j^o j={k~l)i+l 

Using Eq. © one can show that {^^^^\z; N)) = and {^^^\z; N)^^^\z'; N)) = Df5{z - z'), 
where dJJ^ = A^(A^ + 1)D2, D2 = 2A(3s{l - s)T-\ and 6{z) is the Dirac delta function. 
Notice that the first term on the right hand side of Eq. (jSJ is zero due to symmetry. 
Integrating both sides of Eq. over z we obtain 

r]o{z)=exp[2eRX^'\z-N)], (6) 

where x^'^\z; N) = dz' $,^'-'\z' ; N) and we assumed ?7o(0) = 1. According to the cen- 
tral limit theorem x^^\z; N) is a Gaussian random variable with {x^^\z; N)) = and 
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N)) = D^^ z. As a result, the distribution of the sohton amphtude is lognormal 



F(r^o) = (vrI?SV'/'% ' exp [- In^ (r^o) , (7) 

where T)^^ = SD^^^ej^z. The normahzed nth moment of F{riQ) satisfies (r^Q (z))/(?7o(2;))" = 
exp 2n{n — l)D^^\'j^z , from which it follows that the amplitude dynamics is intermittent. 
For a soliton in the zth channel, the amplitude dynamics is given by 

1 dr]i AseRA(3^ , 

+2en^^'\z-N), (8) 

where 



2 \ /T)(0) 



1 ''^^'"'^ 

^^'\z,; N) = -— y J2 sgn(A' - A) 5Z C.'i'- (9) 

iY« i=(fc-i)(j'-j) 



From Eq. © it follows that A^)) = and {C^'\z] N)^^'\z'] N)) = D%^5{z-z'), where 

= [A^(A^ + 1) +i'^]D2- Therefore, the disorder strength is different for different channels. 
To solve Eq. (jH)) we substitute rii{z) = 7]f'\z)7]\^\z) , where Vjf'"' and r/-'^^ represent the drift 
and fluctuating contributions due to the first and second terms on the right hand side of 
the equation, respectively. Assuming that the first term is compensated by appropriately 
adjusting the amplifiers gain, if]f\z) = 1 and r]i{z) = 7]\^\z). As a result, the statistics of 
rji is described by the lognormal distribution (jZj), with V^^j = SD^^ej^z replacing V^^^ 
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The dynamics of the Raman-induced self frequency shift for the reference channel soliton 
is given by 

f3!,'\z) = -^-^ f^z'4{z'). (10) 
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To show that Eq. (jlUp leads to intermittent dynamics for j3^^ we first calculate the nth 
moment: 

X dZi... / dZr,{7]^{Zi)...7]^{Zn)), (H) 

Jo Jo 

where < Zn < ■ ■ ■ < zi < z. Using Eq. © and the fact that the integrals J^' ^ dz' ^^'^\z'; N) 
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are Gaussian random variables that are independent for different i-values we obtain 



X 



where zq 
Kmexp 



n 



m=l 



Thus, (/3o'^'*"(z)) is given by a sum over exponential terms of the form 



8e^ 
/ 15 

2^m — 1 

exp 







32DPel{2m - l)z„ 



(12) 



32m'^D 



(0) 9 



e-hz 



where < m < n and Km are constants. To show intermittency it 
is sufficient to compare the leading term in the sum with the leading term in the expression 
for {I3q\z))^. This calculation yields: 



n\ exp 



32n(n 



l)Dfelz 



Therefore, the leading term in the expression for the normalized nth moment of /Jq**^ is 
exponentially increasing with both z and n^. To illustrate this dynamic behavior we show the 
2;-dependence of the n = 2, 3, 4 normalized moments in Fig. ^for a multichannel transmission 
system with = 50, e/j = 3 x 10""^, s = 1/2, T = 5, and A/5 = 10. These parameters 
correspond to the 101-channel system operating at 10 Gbits/s per channel discussed in 
detail below. One can see that the fourth moment increases much faster with increasing 
z compared with the second and third moments. These results are of special importance 
since other contributions to changes in the soliton parameters are also coupled to the soliton 

is) 

amplitude via integrals over z and thus follow the same statistics as /5q . In particular, 
the s-dependent contribution to the cross frequency shift discussed below and the soliton's 
phase shift, which is given by a^i^z) = 2 dz'rjQ^z'), exhibit similar intermittent dynamics. 
The dynamics of the Raman-induced cross-frequency shift (3q^^ in a two-channel system 



(m 



(13) 



was obtained in Refs. 



15l |. where it was shown that (3^^^ is lognormally distributed. In a 



system with 2N + 1 channels one obtains by a procedure similar to the one used in deriving 
Eq. © 



8eRr]^{zk-i] 



ki 



Z = Zk 



(14) 



Equation (fT^ can be solved by decomposing (3^'^^ into an s-dependent part /3q'^'^^ and an 
s-independent part Pq'^^^ ■ The s-dependent part is given by an integral of 1]^ over z and 

(s) 

thus has similar statistics as (3q . To obtain the evolution of the s- independent part we 



decompose the total disorder field C,''^\z; N) into contributions ^f'\z) coming from different 
channels: ^^^\z; N) = Ei^o sgn(A)d°^(^)> where ^f\z) = ^^Qj/Azi^l Hence, the fields 
satisfy {^f\z)) = and {^l^\z)^l?\z')) = \i\D26ii>6{z - z'). Substituting these relations 
into Eq. (fT^ . going to the continuum limit and integrating over z we obtain 

= -IvH^) E ^ [1 - /'r'(^)] , (15) 



where fii{z) = exp 4sgn{(3i)eRxf\z) , and xf\z) = dz' ^f'\z'). Therefore, f3^^\z) is a 
product of a lognormal variable and a sum over independent lognormal variables. Since the 
terms in the sum on the right hand side of Eq. (fT3j) decrease with frequency as it is 

sufficient to consider only contributions from a few neighboring channels. For a three-channel 
system, for example, Eq. (|THjl simplifies to 

/3f'^(^) = ^[/^-iW-/^i(^)]- (16) 
Using Eq. one finds that the leading term in the expression for the normalized 2nth 



moment of jS^^^^ is exp 32n(n — l)D'^^^e\z /2"', which is exponentially growing with 
Hence, the cross-frequency shift in the three-channel system exhibits intermittent dynamic 
behavior, even though it is not lognormally distributed as in the two-channel case. 

To gain further insight into the intermittent dynamic behavior exhibited by the solitons 
we calculate the nth order equal-distance amplitude correlation functions, which measure 
correlation between amplitudes of solitons from different time slots in the same channel. 
Considering the reference channel we calculate 

Co, = ((r]oo^o,)")/((0(^o"))-l, (17) 

where ?7oo and r^oj stand for the amphtudes of the solitons in the 0th and jth time slots, 
respectively. We still assume that the amplitudes of solitons in other channels are 1. There- 
fore, considering collisions with solitons from the iih channel, for example, the difference 
between the dynamics of the two solitons is due to the fact that the Oj soliton experiences 
the disorder experienced by the 00 soliton with a delay given by /S.Zji = jT/{2f3i). Using 
this fact and the decomposition of into the ^f'\ and assuming j > N, one can show that 
Coj{z) = Voo{z)Voj{^) is lognormally distributed 

exp{-ln'(Co,)/[8elM\S ^ ^ 



)(0).2 



Z. 



8ne%D^N\z) 



Oj 
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where 



bf{z)=AD2 
2D2T{N - i 



K1)|J| 



A/3 



(19) 



for Azji < z < A2;j(i_i). Consequently, the normahzed nth order equal-distance amplitude 
correlation functions are given by 



exp <j 2n^e|j 



D^S\z)-2DPz 



1. 



(20) 



In particular, for z < Azj^ CQj{z) = since the two solitons are uncorrelated. During the 
transient Azjjy < z < Azji the sohtons become correlated due to the effective collision- 
induced disorder. For z > Azji, i.e., after the transient. 



exp 



An^D'^^^elz-8s{l - s^NeH 



-1. 



(21) 



Thus, after the transient the nth order correlation functions grow exponentially with both 
n^ and z, in accordance with the intermittent behavior of the amplitude. Notice that CqJ"* 
decays exponentially with A^|j|, which is the total number of time slots in all other channels 
separating the two solitons. Using Eq. (PT|) with n = 1 one obtains jcor = 1/ (2iVe/j) for the 
typical correlation number. For a 101-channel system operating at lOGbits/s per channel, 
jcor ~ 10^, which is much larger than the number of successive bits that can be corrected 
by current error correction methods (~ 10^), and much smaller than the number of time 
slots that are in transmission in a given channel at any given time (~ 10^). Thus, this type 
of effective collision-induced disorder presents a challenge for conventional error correction 
methods. 

We now relax the frozen disorder assumption and take into account the dynamics of 
soliton amplitudes in all channels. In this case Eq. (jlj) is replaced by 

1 dr]o isenAp 



rjo dz T 
+2e^e(°)(^;iV), 

where rii{z) = exp [2eji dz' ^^^\z'] N)'j and 

i^'\z;N) 



\i\r]i{z)sgn{(3i 



(22) 



-AiYy]^*Wsgn(A 

Az^c' 
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) E 

j=(fc-i)i+i 



(23) 
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Expressing rio{z) as a product of an s-dependent and an s-independent parts: rio{z) 



{z)r]Q ' {z) ^ substituting into Eq. 



and integrating over z we obtain 



id)( \ 

Vo (^) = exp 



T 



dz'r]i{z') 



(24) 



and 



r^(^)(z) 



exp 



2en / d^'eW(z';iV) 



(25) 



It follows that r]^\z) is no longer deterministic. Moreover, since r]^\z) is proportional 
to the exponent of the integral over //^(-z), where r]i{z) is lognormal, one can expect the 
departure of the t]^'' statistics from Gaussian statistics to be stronger than lognormal. Con- 
sider now the statistics of 'r]"Q \ Using Eq. (j23|) one can show that {i^^\z]N)) = and 
(^^(2; N)i^''\z'; N)) = Df{z)5{z - z'), where 



i\ exp 



SD^^elz 



(26) 



and we assumed {rii{z)riii{z')C,ijC,iiji) = {rii{z)rii'{z')){(ij(i'ji). As a result, the T^g'^'^-distribution 
is the lognormal distribution given by Eq. ((Zj) with V^^^ replaced by V^^\z), where 



exp[8D'i^elz 



(27) 



A direct consequence of Eq. ()27|) is that the nth moments of the 77Q -distribution are super- 
exponentially increasing with z, although the factors d>D'^^e\z are much smaller than 1. 

From the practical point of view it is important to understand the influence of the inter- 
mittent dynamic behavior of the soliton parameters on the BER. The contribution of the 
collision-induced pulse decay to the BER was discussed in detail in previous works (see Ref. 
jlT^ and references therein). Moreover, the small-77 tail of the lognormal distribution lies 
below the corresponding tail of the Gaussian distribution, whereas the large-77 lognormal 
tail lies above the corresponding Gaussian tail. As a result, strong effects due to deviations 
from Gaussian statistics are related to relatively large ?7-values. When the position dynam- 
ics or the frequency dynamics are coupled to the amplitude dynamics such large r^-values 
can lead to significant increase in the BER due to walk-off of the soliton from its assigned 
time slot. Therefore, we focus our attention on contributions to the BER due to the large-77 
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lognormal tail. We consider a 101-channel system operating at lOGbits/s per channel and 
emphasize that state-of-the-art experiments with dispersion- managed solitons demonstrated 
multichannel transmission with 109 channels at 10 Gbits/s per channel over a distance of 
2 X 10^ km 181. We use the following parameters, which are similar to the ones used in 

n 

multichannel soliton transmission experiments [19[. Assuming that T = 5, A/5 = 10 and 
s = 1/2, the pulse width is 20 ps, eR = 3x 10~^, the channel spacing is 75 GHz, and D2 = 1. 
Taking j32 = — Ips^/km, the soliton-peak-power is Pq = 1.25 mW. For these values the dis- 
order strength is I)^q{z) = 1.8 x 10~^z for the reference channel and V^^q\z) = 3.6 x 10~^z 
for the two outermost channels. For z = 25, corresponding to transmission over 2 x 10^ km, 
^(^(25) = 0.046 and V^^o\25) = 0.091. 

For this system we evaluated the contributions to BER due to the Raman-induced cross 
frequency shift, the Raman-induced self frequency shift and the "ideal" component of the 
collision-induced position shift, i.e., the position shift due to soliton collisions in the absence 
of perturbations. The calculations show that the dominant contribution to the BER is due 
to the s-dependent part of the cross- frequency shift /3q'^'^'*. The position shift induced by (3^^'^^ 
is obtained by taking the continuum limit in Eq. ()14j) and integrating the s-dependent term 
twice with respect to z: 

yr(^)=-^rd.'r'dA„^(."). (28) 



^0 

~(cd 



3T 

The position shift with a fixed amplitude rio{z) = 1 is yl)'"''{z) = — [32 Nenszj) / {3T) and 
the relative position shift is AyQ^'^\z) = y^'^\z) — y^Q'^\z). We assume that y^^'^^ can be 
compensated by employing filters. Therefore, the total energy of the soliton at a distance z 
is 

1(7)0, Ayf = r)l / dt cosh-2[r/o(t - Ayi^"^)]. (29) 

J-T/2 

Occupied time slots are considered to be in error, if J(?7o, A?/q'^'^'') < I{z = 0)/2 ~ 1. To 
estimate the BER we numerically integrate Eq. ()28|) coupled to Eq. ® for different realiza- 
tions of the disorder $,^^\z, N) and calculate the fraction of errored occupied time slots. The 
BER in a generic channel i 7^ is calculated in a similar manner, where r]Q and ^^^\z,N) 
are replaced by rji and ^^'^\z, N), respectively. Figure |21 shows the 2;-dependence of the BER 
in channels i = (the reference channel), i = 25, and i = 50 (the outermost channel) for 
the aforementioned system. One can see that the BER in the reference channel increases 
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from values smaller than 10~^ for z < 15 {x < 1.2x 10^ km) to about 8.2 x 10~^ at z = 25.0 
(x = 2 X 10^ km). Furthermore, for intermediate distances 15 < z < 20, the BER value in 
the outermost channels can exceed that in the reference channel by several orders of magni- 
tude, even though, the disorder strengths differ by only a factor of 2. This behavior presents 
another challenge to conventional error correction methods based on knowledge gained from 
single- or few-channel transmission systems. 

To better understand error generation due to (3^^'^^ we analyzed the 2;-dependence of 
contributions to the BER coming from different regions in the rjQ — Ayl^'^^ plane. The results 
of this analysis are presented in Fig. (jH}. At 2 = 15 the dominant contribution to the BER 
comes from the domain rjo < 0.7 and ^u^q'^^ > 0, i.e., from decaying solitons with relatively 
large positive values of Ay^^^K For z > 17 the dominant contribution comes from the region 
rjo > 1.0 and Ayl^'^^ < 0, which corresponds to solitons with relatively large amplitudes 
and large negative values of Ay^'^\ The latter contribution is associated with the large- 
Tj lognormal tail of the amplitude distribution. Figure IH shows the mutual distribution 
function G{r]Q, Ay^^^'^^) at 2; = 25 and the two domains giving the main contributions to 
the BER. It can be seen that this distribution is very different from the one observed for 
single-channel soliton propagation in the presence of amplifier noise (see Fig. 1 in Ref. P|). 
While the latter distribution is approximately symmetric about Ayo = and 770 = 1, the 
former is strongly asymmetric with an extended tail in the large-770 and large-negative-A?/Q'^'^^ 
region. The strong asymmetric form of G{i]q, AyQ^^^) in our case is due to the strong coupling 
between the position dynamics and the amplitude dynamics and the lognormal statistics of 
the soliton amplitude. Thus, we find that amplitude dynamics plays a dominant role in error 
generation in massive multichannel optical fiber transmission systems, a situation which is 
very different from the one observed in single-channel transmission systems j3| . 

In summary, we studied soliton propagation in massive multichannel optical fiber commu- 
nication systems taking into account the effects of delayed Raman response and the random 
character of pulse sequences. We found that the soliton parameters exhibit intermittent 
dynamic behavior and showed that the cross frequency shift is the main mechanism leading 
to bit pattern deterioration and to relatively large values of the bit-error-rate. We empha- 
size that similar dynamic behavior is expected in massive dispersion-managed multichannel 
transmission systems as well. In such systems the Raman-induced energy exchange in col- 
lisions will lead to lognormal statistics for the pulse amplitudes. In addition, the frequency 
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and position dynamics will be affected by a variety of amplitude-dependent perturbations 
due to Kerr nonlinearity. The coupling of the frequency and position dynamics to the am- 
phtude dynamics will lead to intermittent dynamics of the pulse frequency and position and 
to relatively large values of the bit-error-rate. 



[1] G. P. Agrawal, Nonlinear Fiber Optics (Academic, San Diego, CA, 2001). 
[2] C. R. Menyuk, Opt. Lett. 20, 285 (1995). 

[3] G. Falkovich, I. Kolokolov, V. Lebedev, V. Mezentsev, and S. Turitsyn, Physica D 195, 1 
(2004). 

[4] S. A. Derevyanko, S. K. Turitsyn, and D. A. Yakushev, J. Opt. Soc. Am. B 22, 743 (2005). 
[5] E. T. Spiller, W. L. Kath, R. O. Moore, and C. J. McKinstrie, IEEE Photon. Technol. Lett. 
17, 1022 (2005). 

[6] U. Frisch, Turbulence: The Legacy of A. N. Kolmogorov, (Cambridge University Press, Cam- 
bridge, England, 1995). 
[7] F. M. Mitschke and L. F. Mollenauer, Opt. Lett. 11, 659 (1986). 
[8] J. P. Gordon, Opt. Lett. 11, 662 (1986). 
[9] S. Chi and S. Wen, Opt. Lett. 14, 1216 (1989). 
[10] B. A. Malomed, Phys. Rev. A 44, 1412 (1991). 

[11] C. Headley III and G. P. Agrawal, J. Opt. Soc. Am. B 13, 2170 (1996). 

[12] S. Kumar, Opt. Lett. 23, 1450 (1998). 

[13] T. L. Lakoba and D. J. Kaup, Opt. Lett. 24, 808 (1999). 

[14] A. Peleg, Opt. Lett. 29, 1980 (2004). 

[15] Y. Chung and A. Peleg, Nonlinearity 18, 1555 (2005). 

[16] The dimensionless z in Eq. is z = (/32x)/(2to), where x is the actual position, tq is the 
sohton width, and P2 is the second order dispersion coefficient. The dimensionless retarded 
time is t = t/tq, where r is the retarded time. The spectral width is uq = 1/{it'^tq) and the 
channel spacing is Av = {t:/S.(3uq) /2. ^ = Ej^JP^^^ where E is the electric field and Pq is the 
peak power. The dimensionless second order dispersion coefficient is d = — 1 = /?2/(7-Po''"o )> 
where 7 is the Kerr nonlinearity coefficient. The dimensionless Raman coefficient is e^j = 
0.006/ro, where tq is in picoseconds. 



12 



[17] F. Forghieri, R. W. Tkach, and A. R. Chraplyvy, in Optical Fiber Telecommunications III, 
edited by 1. P. Kaminow and T. L. Koch (Academic, San Diego, CA, 1997), Chap. 8, Sec. 
VIII. 

[18] L. F. Mollenauer, A. Grant, X. Liu, X. Wei, C. Xic, and I. Kang, Opt. Lett. 28, 2043 (2003). 
[19] L. F. Mollenauer and R V. Mamyshev, IEEE J. Quantum Electron. 34, 2089 (1998). 



13 




FIG. 1: Normalized moments of the reference-channel soliton self frequency shift (/3q )/(/3o )" '^s 
propagation distance z for a 101-channel system operating at 10 Gbits/s per channel. The solid, 
dashed, and dotted lines correspond to the n = 2, 3, 4 normalized moments calculated by using Eq. 
((T^. respectively. 
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FIG. 2: BER vs propagation distance z for a 101-channel transmission system operating at 10 
Gbits/s per channel. The squares, circles, and triangles represent the BER at channels z = 
(central), i = 25, and i = 50 (outermost), respectively. 
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FIG. 3: The z-dependence of different contributions to the BER for the reference channel in a 
101-channel system operating at 10 Gbits/s per channel. The up triangles correspond to the total 
BER. The squares, circles, and down triangles correspond to contributions coming from the regions 
rj > 1.0 and y <0, r] < 0.7 and y > 0, and 0.7 < r] <1.0, respectively, in the r/o — ^y^^ plane. 
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FIG. 4: Mutual distribution function G{r)o,AyQ for a reference channel soliton at z=25 in a 
101-channel system operating at 10 Gbits/s per channel. The arrows A and B show the domains 
giving the main contributions to the BER. 
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